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Abstract —In the current DVB generation, satellite terminals 
are expected to be interactive and capable of transmission in the 
return channel with satisfying quality. Considering the bursty 
nature of their traffic and the long propagation delay, the use of a 
Random Access technique is a viable solution for such a Medium 
Access Control (MAC) scenario. In this paper Random Access 
communication design in DVB-RCS2 is considered with partic¬ 
ular regard to the recently introduced Contention Resolution 
Diversity Slotted Aloha (CRDSA) technique. The paper presents 
a model for design and tackles some issues on performance 
evaluation of the system by giving intuitive and effective tools. 
Moreover, dynamic control procedures able to avoid congestion at 
the gateway are introduced. Results show the advantages brought 
by CRDSA to DVB-RCS2 with regard to the previous state of 
the art. 

Index Terms —Retransmission, Successive Interference Cancel¬ 
lation, Slotted Aloha, Random Access, Congestion Avoidance, 
Dynamic Control, Satellite Broadcasting, Interactive Terminals 


I. Introduction 

Current DVB standards such as Digital Video Broadcasting 
over Satellite are increasingly going towards the need of 
high-speed bidirectional networks for consumer interactivity 
|Tj. Nevertheless in consumer type of interactive satellite 
terminals, users generate a large amount of low duty cycle 
and bursty traffic with frequent periods of inactivity in the 
return link. Under these operating conditions, the traditionally 
used Demand Assignment Multiple Access (DAMA) satellite 
protocol does not perform optimally, since the response time 
for the transmission of short bursts can be too long. For this 
reason, in the recently approved specification for the next 
generation of Interactive Satellite Systems (DVB-RCS2) [2], 
the possibility of sending logon, control and even user traffic 
using Random Access (RA) in timeslots specified by the 
Network Control Center (NCC) is provided to Return Channel 
Satellite Terminals (RCST). In particular, two methods are 
considered for RA: the first one is Slotted Aloha (SA) (3) 0, 
the second is called Contention Resolution Diversity Slotted 
Aloha (CRDSA) 0. 

SA represents a well established RA technique for satellite 
networks in which users send their bursts within slots in a 
distributed manner, i.e. without any central entity coordinating 
transmission. This allows to reach an average throughput 
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around 0.36 [packets/slot] despite the possibility of collision 
among bursts from different users. Nevertheless, in prac¬ 
tice SA works with very moderate channel load to ensure 
acceptable delay and loss probability (PER = 10“ 3 for 

10 -3 [packets / slot]). 

This gap has been recently filled with the introduction of 
a new technique named CRDSA that is able to boost the 
throughput even up to values close to 1 [packets/slot]. This 
technique is based on the transmission of a chosen number 
of replicas for each burst payload as in Diversity Slotted 
Aloha (DSA) (6j. However, differently from DSA, in CRDSA 
each burst copy has a pointer to the location of the other 
replicas of the same burst payload so that when the frame 
arrives at the receiver a Successive Interference Cancellation 
(SIC) process is accomplished. The SIC process consists in 
removing the content of already decoded bursts from the slots 
in which collision occurred. Doing so, it is possible to restore 
the content of bursts that had all their replicas colliding thus 
boosting the throughput. Nevertheless the use of interference 
cancellation can benefit from power unbalance among different 
sources 0 ® 

Considering these recent findings on Aloha-based RA and 
the fact that in © contention control mechanisms to avoid 
congestion are claimed to be out of scope for the document, 
in this paper a complete model for the design of a congestion- 
free systemusing CRDSA is presented. To do so, we rely on 
the analysis we carried out in (9), in which the need for such 
a model was firstly introduced. 

The main contributions of this paper are: the presentation 
of a stability model and of some crucial evaluation parameters 
(packet delay and FET) developed in [TO]] for SA and here 
adapted to CRDSA; a thorough analysis of the design aspects 
to take into account in such a system; the adaptation of 
the dynamic control limit policies developed for SA and the 
analysis of their convenience when compared to static CRDSA 
design; the extension of all these results to the case of infinite 
population that is usually not considered even though in some 
scenarios is the most appropriate population representation. 

The remainder of this paper is organized as follows. In 
Section II a brief overview of the CRDSA scheme is given and 
the congestion problem that arises when introducing retrans¬ 
missions in such a scheme is stated. Section III presents the 
stability model used for analysis and design in the subsequent 
sections and defines the concept of stability. Section IV regards 
the definition of a Markov chain for packet delay computation 
and its validation. In Section V the definition of First Exit Time 
is given and the formulas for its calculation are presented; 
moreover the adaptation of the model for computation reduc¬ 
tion and for the extension to the case of infinite population are 
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discussed. Section VI exploits the model and the evaluation 
parameters introduced in the previous sections by giving a 
general overview of how they are used in the design phase. 
Section VII compares the advantages and disadvantages of 
various Aloha-based RA techniques and packet degree choices. 
Finally Section VIII introduces some dynamic control policies 
already used in m for SA while Section IX and X discuss 
their application respectively to the case of finite and infinite 
population. Section XI concludes the paper. 

II. CRDSA System Overview 

Consider a multi-access channel populated by a total num¬ 
ber of users M. Users are synchronized so that the channel is 
divided into slots and N s consecutive slots are grouped in a 
frame. The probability that at the beginning of a frame an idle 
user will send a new packet is po. When a frame starts, users 
willing to transmit place d copies of the same packet over the 
N s slots of that frame. The number of copies d can be either 
the same for each packet or not. The first case is known as 
Constant Replication CRDSA (CR-CRDSA); the second case 
is known as Variable Replication CRDSA (VR-CRDSA). In 
this second case the number of copies is randomly generated 
using a pre-determined probability mass function ED 

Packet copies are nothing else that redundant replicas except 
for the fact that each one contains a pointer to the location 
of the others. These pointers are used in order to attempt 
restoring collided packets at the receiver by means of SIC, 
i.e. by removing the interfering content of already decoded 
packets thanks to their location’s knowledge. Nevertheless, 
the physical structure of the burst copies used is equal to 
the one used for dedicated access with respect to the burst 
construction, the waveform, the timeslot structure and the burst 
reception. 

Throughout the paper perfect SIC and ideal channel are 
assumed, which means that the receiver is always able to 
perfectly apply SIC without errors and no external sources of 
disturbance such as noise that might require SNR estimation 
ED are present. Therefore, the only cause of unsuccessful 
decoding of packets is interference among their[] Given these 
conditions, consider Figure [T] representing a practical example 
of SIC. Each slot can be in one of three states: 

• no packet’s copies have been placed in a given slot, thus 
the slot is empty; 

• only 1 packet’s copy has been placed in a given slot, thus 
the packet is correctly decoded; 

• more than 1 packet’s copy has been placed in a given 
slot, thus resulting in collision. 

In DSA, if all copies of a given packet collided the packet is 
surely lost. In CRDSA, if at least one copy of a certain packet 

! For the case in which non-ideal channel estimation is considered, the 
reader can refer to [5] where it has been demonstrated that there is no 
appreciable impact on the CRDSA performance for moderate SNR. Moreover, 
as suggested by one of the reviewers that we kindly thank, we specify here 
that also situations such as outage that affect the number of backlogged users 
are not be a problem as long as the channel is nominally stable. In fact, if 
there is only one point of stability, the channel will rapidly set back to the 
operational point after the outage period. This sentence is anticipated here but 
will be clearer in the subsequent sections. 
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Figure 1. Example of frame at the receiver for CRDSA with 2 copies per 
packet; plain slots indicate that a transmission occurred for that user in that 
slot. 

has been correctly received (see User 4), the contribution of 
the other copies of the same packet can be removed from the 
other slots thanks to the knowledge of their location provided 
by the pointers contained in the decoded packet. This process 
might allow to restore the content of packets that had all their 
copies colliding (see User 2) and iteratively other packets may 
be correctly decoded up to a point in which no more packets 
can be restored (see User 1 and 3) or until the maximum 
number of iterations Imax for the SIC process is reached, if a 
limit has been fixed. 

If a packet has not been successfully decoded at the end of 
the SIC process and a retransmission request is sent by the 
NCC, a feedback is created between the transmitter and the 
receiver side. This feedback could cause the channel to get to a 
state of congestion in which the load is too big because of too 
many packets to retransmit and causes excessive packet loss 
due to simultaneous destructive transmissions. For this reason, 
a careful design of the system based on a stability model is 
of great interest for practical DVB-RCS applications. 

III. Stability Model 

Consider the RA communication system presented in Sec¬ 
tion [II] Each user can be in one of two states: Thinking (T) 
or Backlogged (B). Users in T state generate a new packet for 
transmission with probability po over a frame interval; if so, 
no other packets are generated until successful transmission 
for that packet has been acknowledged. Users in B state have 
unsuccessfully transmitted their packet and keep attempting to 
retransmit it with probability p r over a frame interval. In the 
followings, we assume that users are acknowledged about the 
success of their transmission at the end of the frame in which 
the packet has been transmitted (i.e. immediate feedback). 
Nevertheless this constraint will be relaxed in the last sections. 

A. Equilibrium Contour 

Defining 

• Ng : backlogged users at the end of frame / 

• Gg = Nb n Vr '• expected channel loacj^Jof frame / due 
to users in B state 

2 In this paper, with channel load we refer to the logical channel load that 
is the number of packets sent normalized over the frame size, regardless of 
the number of copies sent per packet. 
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• Gtp : expected channel load of frame / due to users in 
T state 

• Gj N = Gj, + Gg : expected total channel load of frame 

/ 

• PLRf ( Gj N , N s , d, Imax) • expected packet loss ratio of 
frame / 

• Gqut = GjnO- — PLR(G{ N ,N s ,dJ max )) : part of 
load successfully transmitted in frame /, i.e. throughput. 

The equilibrium contour can be written as 


Gt = Gout = G/at(1 — PLR(Gin , N s , d, Imax)) (1) 

that represents the locus of points for which at any time the 
expected channel load due to users in T state is equal to the 
expected throughput Gout • Notice that the frame number / 
has been omitted, since in equilibrium state this condition is 
expected to hold for any frame. 

From the definition above we can also gather that the 
expected number of backlogged users remains the same frame 
after frame. Therefore 

AT _G in PLR(G in ,N s , d, Imax)Us 

=- U) 

Pr 

Equations 0 and © completely describe the equilibrium 
contour on the ( GtMb) plane for different values of G/a^] 


B. Channel load line 

Notice that the equilibrium contour itself is an infinite set 
of equilibrium points. In fact, we have not considered M and 
Po so far. This two essential parameters constitute the channel 
load line and provide the knowledge of the actual stability 
points for a certain scenario. 

Consider M and po to be constant (i.e. stationary input). 
Given a certain number of backlogged packets, the channel 
load line expresses the expected value of channel load due to 
users in T state. For the finite population case, the channel 
load line can be defined as 


Gt 


M -N b 
N s 


Po 


(3) 


while for M -A oc the channel input can be described as a 
Poisson process with expected value A [thinking users] 0 so 


that 



(4) 


for any N B (i.e. the expected channel input is constant and 
independent of the number of backlogged packets). 


3 Concerning the values used for the Packet Loss Ratio, it is known from 
the literature (5j that the relation between PLR(Gin) and Gjn can not be 
tightly described in an analytical manner. For this reason PLR values used 
in this work are taken from simulations. However, once simulations for the 
open loop case are accomplished, they can be used in the stability model to 
see how the performance changes when modifying some crucial parameters 
such as p r , po and M without the need of running brand new simulations. 
Moreover simulated values can be substituted with tight analytical equations 
in case any will be found in the future. 



Po = 0.143 
p r = 0.5 

M = 350 


(Gr,Ng) = (0.49,8.2) 



(b) Unstable channel (finite M) 

Po = 0.143 
Pr = 1 

M = 350 

(Gf 1 ,^! 1 ) = (0.495,3.9) 
(G%,NU) = (0.44,42.7) 
(Gp,N§ 2 ) = (6 • 10 _3 ,346) 



(c) Unstable channel (infinite M) 

A = 40 
p r = 0.5 

M oo 

(Gf. 1 ,-/Vg 1 ) = (0.4,1.7) 
(Gj,,Ng) = (0.4,107.32) 
(Gp,N§ 2 ) = (0.4, oo) 



(d) Overloaded channel 
Po = 0.18 
Pr = 1 

M = 350 

(Gf,V|) = (6 • 10 -3 ,346) 


Figure 2. Examples of stable and unstable channels for CRDSA (2 replicas) 
with N s = 100 and Imax = 20. Stable equilibrium points are marked with 
a black dot. 


C. Definition of stability 

Figure [2] shows the equilibrium contour and the channel 
load line for various cases. The equilibrium contour divides 
the ( N b ,Gt) plane in two parts and each channel load line 
can have one or more intersections with the equilibrium 
contour. These intersections are referred to as equilibrium 
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points since Gout = Gt- The rest of the points of the 
channel load line belong to one of two sets: those on the 
left of the equilibrium contour represent points for which 
Gout > Gt, thus situations that yield to decrease of the 
backlogged population; those on the right represent points for 
which Gout < Gt , thus situations that yield to growth of 
the backlogged population. 

Therefore we can gather that an intersection point where 
the channel load line enters the left part for increasing 
backlogged population corresponds to a stable equilibrium 
point. In particular, if the intersection is the only one, the 
point is a globally stable equilibrium point (indicated as 
G^Me) while if more than one intersection is present, it 
is a locally stable equilibrium point (indicated as G^,7Vj). 
If an intersection point enters the right part for increasing 
backlogged population, it is said to be an unstable equilibrium 
point (indicated as G^,A/^). 


D. Channel state definition 


Figure |2(a)| shows a stable channel. The globally stable 
equilibrium point can be referred as channel operating point 
in the sense that we expect the channel to operate around that 
point. With the word around we mean that due to statistical 
fluctuations, the actual Gt and Gb (and thus also Gin and 
Gout ) may differ from the expected values. In fact, the actual 
values have binomially distributed probability for Gb and 
Gt with finite M and Poisson distributed probability for Gt 
with infinite M. Nevertheless, the assumption of taking the 
expected value has been already validated in j9j. 

Figures |2(b) and 2(c)| show two unstable channels respec¬ 
tively for finite and infinite number of users. Analyzing this 
two figures for increasing number of backlogged packets, the 
first equilibrium point is a stable equilibrium point. Therefore 
the communication will tend to keep around it as for the 
stable equilibrium point in Figure 2(a) and we can refer 
to it once again as channel operating point. However, due 
to the abovementioned statistical fluctuations, the number of 
backlogged users could pass the second intersection and start 
to monotonically increase. 


In the case of finite M, this increment goes on till a new 
intersection point is reached, while in the case of infinite M 
the expected number of backlogged users increases without 
any bound. In the former case, this third intersection point is 
another stable equilibrium point known as channel saturation 
point , so called because it is a condition in which almost any 
user is in B state and Gout approaches zero. In the latter 
case, we can say that a channel saturation point is present for 

N b oo. 

Finally Figure [2(d)] shows the case of an overloaded channel. 
In this case there is only one equilibrium point corresponding 
to the channel saturation point. Therefore, even though the 
channel is nominally stable, the point of stability occurs in 
a non-desired region. For this reason this case is separated 
and distinguished from what is intended in this work as stable 
channel. 


IV. Packet Delay Model 

Once we are sure that the channel is working at its channel 
operating point , we would like to know what is the expected 
distribution and average delay associated to packets that are 
successfully received at the gateway. This can be described 
using a discrete-time Markov chain with the two states B and 
T (Figure [3]). 

The edges emanating from the states represent the state tran¬ 
sitions occurring to users which depend on p r and the expected 
PLR at the channel operating point. If a packet transmitting 
for the first time receives a positive acknowledgment the user 
stays in T state while in case of negative acknowledgment 
the user switches to B state until successful retransmission 
has been accomplished. We assume a frame duration to be 
the discrete time unit of this Markov chain. Therefore the 
packet delay D p kt is calculated as the number of frames that 
elapse from the beginning of the frame in which the packet 
was transmitted for the first time, till the end of the one in 
which the packet was correctly received. 



According to the definition given above, the delay distribu¬ 
tion is entirely described by 


Pr{D pkt = /} 


1 — PLR, for / = 1 

PLR [p r (1 — PLR)]- 

•[1 — Pr + PLR PrY~ 2 , for / > 1 

(5) 


which results in the average delay 


Av[D pkt ) = ]T/- Pr{D pkt = /} (6) 

f= 1 

Equation [6] can alternatively be rewritten in a simpler and 
more intuitive form by means of Little’s Theorem |T0| GD 
considering that the average number of backlogged users in a 
stable system is equal to the average time spent in backlogged 
state, multiplied by the arrival rate of new packets Gt (that we 
know to be equal to Gout at the operational point). Therefore 

Av[Dpkt] = gTYUNs (7) 

where the presence of N s in the formula has the aim of 
normalizing the delay to the frame unit. 
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Pr 

Av \Dpkt\ana 

Av [Dpkt]sim 

PLR 

1 

543.48 

541.89 

0.99816 

0.4 

3.17 

3.16 

0.465 

0.2 

2.01 

2.12 

0.168 

0.05 

2.53 

2.66 

0.0713 

0.005 

5.08 

5.22 

0.02 


Table I 

Analytic and simulated average packet delay at the channel 
OPERATING POINT FOR CRDSA WITH N s = 100 slots, Imax = 20, 
M = 350 AND po = 0.18. 


A. Model validation and retransmission probability's effect 

Figures]?] and [5] show some results based on the example of 
overloaded channel given in Section III, in order to validate 
the model and demonstrate the effect of p r on the packet delay. 
In particular p r is progressively decreased in order to evaluate 
how distribution and average packet delay change. Figure [4] 
represents the channel load line and equilibrium curves for 
several p r values. 

When p r is decreased, the equilibrium contour moves up¬ 
wards and the point of equilibrium is found for bigger values 
of throughput while N B decreases. Therefore, from Little’s 
theorem, it is intuitive to expect that users will spend less 
time in B state. If we keep decreasing p r after the throughput 
peak is reached, the throughput starts to decrease while the 
number of backlogged users increases again thus yielding to 
an increase of the average packet delay. This is confirmed by 
the results given in Table [I] 

Figure [5] shows that while having a small p r is beneficial for 
the PLR , it also means that users in B state will generally wait 
longer before retransmitting a packet. As a consequence, the 
smaller the p r value the less steep the cumulative distribution 
curve is, so that at a certain point curves with bigger p r 
overcome it in terms of probability to receive a packet before 
a certain time deadline. Finally the simulated average packet 
delay associated to the examples in Figure [5] is reported in 
Table [I] together with the expected results, demonstrating the 
validity of the model. 

V. First Exit Time 

In Section III the concept of unstable channel has been 
introduced. Beside knowing that a certain channel is unstable, 
it is also of big interest to quantify how much the channel is 
unstable. The First Exit Time (FET) calculates the time spent 
around the operational point by an unstable channel before 
leaving the stability region. To calculate this parameter we 
consider a discrete-time Markov chain which uses the number 
of backlogged users N B as describing state, since equilibrium 
points and more generally the state of the communication can 
be uniquely identified with it. 

A. Transition matrix definition 

Let us define P as the transition matrix for this Markov 
chain, with each element pij denoting the probability to pass 
from a number of backlogged users N B — j to a number of 



Figure 4. Equilibrium curves for CRDSA with N s = 100 slots, Imax = 
20, M = 350 and p 0 = 0.18. 



Figure 5. Packet delay cumulative distribution for CRDSA with N s = 
100 slots, Imax = 20, M = 350 and po = 0.18. 


backlogged users N B = i. To calculate each p^ value we 
need to know the probability that having a certain number of 
backlogged users j and a certain number of thinking users 
M — j, respectively b backlogged users and t thinking users 
have transmitted and s = t — (i — j) out of t + b packets have 
been correctly received. Therefore 

t =o b =o v J (8) 

• {fyp b r0 ~ PrY~ b ‘ q(s\t + b) 

where q(s\t-\-b) is the probability that s packets are correctly 
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Figure 6. P matrix for CRDSA with N s = 100, Imax = 20, M = 350, 
po = 0.143, p r = 1. 

decoded if £ + b packets are transmitted Let us then call Q 
the matrix with elements corresponding to probability 

q(s\t + b). 

In Figure [6] an example of the resulting P matrix is shown. 
Accordingly to the bidimensional representation used in Sec¬ 
tion III, we can see here that there are two peaks corresponding 
to the two expected locally stable equilibrium points. However, 
it is more fruitful for our considerations to look at this graph 
from an (i, j) plane perspective plotting a line for i — j. 
Points below this line represent decrements of the number of 
backlogged packets while points above represent increments. 
Non-zero probability points look like a wake close to the line 
i = j, highlighting that possible transitions are found only for 
small changes in the number of backlogged packets compared 
to the total population M. Taking single values of j, we can 
roughly identify 3 regions of the wake: 

• for small values of j (close to the channel operating point) 
non-zero probabilities are found both for increase and 
decrease of Nb\ 

• for values of j close to M (i.e. close to the channel 
saturation point) non-zero probabilities are found both 
for increase and decrease of Nb’, 

• for the rest of the j values non-zero probabilities are 
almost exclusively found for i > j so that the number 
of backlogged users monotonically increases. 

Therefore in this particular case the system has probability 
zero of returning to the stability region after Nb crosses 120. 
We want to highlight here that the presence of a no-return 
point is scenario-dependent in the sense that for a different 
set of values M, p 0 » Pr , N s it is not necessarily true that 
the probability of coming back from a certain value N B is 
zero: sometimes could be negligible, sometimes could not, 

4 For the same reason claimed in footnote 3 for the PLR, also q(s\t + b) 
values are calculated from simulations over a big amount of runs (over 10 6 per 
(t + b) value) since analytical formulas able to tightly describe the probability 
for a packet to be correctly decoded have not been found yet. Nevertheless 
in the followings a method to reduce this computational effort is introduced. 



Figure 7. Evolution of B? for the P matrix in figure [6] 

even though for the typical settings used in these scenarios the 
last case is really rare. Nevertheless this highlights one of the 
advantages of calculating the P matrix: it tells us theoretically 
if there is any probability of coming back to the region of 
stability and in the case the answer is no, tells us what is the 
no-return value of Nb- 

B. State probability matrix and FET calculation 

Based on the P matrix, it is now possible to calculate the 
probability P/ +1 to be in a certain state Nb = i at frame /+1, 
given a certain probability distribution B? = [Bq B{ • • • B^ M ] 
at frame /: 


■r/+ 1_ 


~Po,o 

P0,M " 


AT 

b( +1 

= 

Pl,0 



B{ 

d/+! 

U M 


Pm, o * 

• * Pm,m 




Generalizing Equation [9] the state distribution at each frame 
(/ + 1) can always be calculated from the initial probability 
distribution B° = [#§ T>o Bq 4 } an d the matrix P^ +1 . 
Given the state probability matrix, it is now possible to 
calculate the average and cumulative distribution of the FET 
from a given initial state. In this example we consider as initial 
state B° = [ 1 0 0 • • • 0 ] that is the case in which the 
communication starts with zero backlogged users. 

Figure [7] reports the evolution frame after frame of the state 
probability for the P illustrated in Figure [6] As we can see, 
when / increases the probability of being around the locally 
stable point in the region of high throughput decreases while 
the probability of being in saturation increases. 

C. Reduced transition matrix 

The entire description of any state of this Markov chain is 
of great interest but most of the times impractical, mainly 
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Figure 8. FET cumulative distribution for CRDSA with N s = 100, I max = 
20, M = 350, po = 0.18, p r = 1. 



Figure 9. Simulated and analytic average FET when increasing N^ bs for 
CRDSA with N s = 100, Imax = 20, M = 350, po = 0.18, p r = 1. 


because the dimension of the P matrix is equal to M 2 . 
Therefore it could take a huge computational effort to calculate 
it. Nevertheless it is of interest to know the probability that 
the stability region around the operational point has been left 
in a certain moment, regardless of the actual state once we 
enter the region that yields to saturation. 

In fact, it has been already pointed out in the description 
of Figure [6] that the further the communication gets from the 
first point of local stability, the smaller is the probability that 
the communication will come back to the stability region. For 
this reason, as a first approximation, it is possible to consider 
an absorbing state N^ bs = \Ng~\ + 1 that groups all the 
probabilities for states from \Ng~\ + 1 to M. Doing so, our 
P matrix can be rewritten as 


P = 


Pop 


Pn%, o 


1 - (Ea=0P«,0) 


Po,N% 


Pl,N% 


1 - (J2a=0Pa,Nu) 


( 10 ) 

The P matrix in Equation [lO] takes much less effort than 
calculating the entire P matrix since only (\Ng~\ + l) 2 
elements need to be computed thus speeding up calculations. 

_ N U 

Considering that J2a= o P<*,j must be equal to 1, the probability 
to get to the absorbing state Ng bs from a certain state j can be 

computed as pij = 1 — (J2a=oPa,j)- Moreover, since Ng bs 
is an absorbing state, any probability to leave this state is zero 
while the probability of staying in the absorbing state is equal 
to 1. 


validity of this approximation. Figure [8] shows the distribution 
for the results obtained with the reduced matrix and the 
simulated results. Regarding simulations, if N B crosses the 
value \Ng ] but it comes back to the stability region, the 
FET is reset and the next exit is waited to set it again. For 
the analytic results the FET is calculated accordingly to the 
reduced matrix in Equation [TO] 

As already pointed out in the description of Figure [6] 
when the number of backlogged packets N B is still close to 
the unstable equilibrium point, a certain probability that the 
communication will come back from the instability region is 
present. This is the reason why in Figure [8] the two curves 
show a noticeable difference. However, this problem can be 
solved by considering N^ bs = \N B ] + 1 + A as absorbing 
state, where A is a positive integer big enough to ensure that 
the probability of false exits from the stability region (i.e. the 
probability that N B crosses the stability region but comes back 
to values < N^ bs ) is sufficiently low. Doing so, simulation and 
analytical results with the approximation of absorbing state 
perfectly match. Figure [9] shows how the difference among 
the two average FETs decreases while increasing the chosen 
value for N^ bs . 

In conclusion, we can state that analytic results with the 
approximation of absorbing state represent a lower bound 
for the actual FET. Therefore during this kind of analysis 
a tradeoff between tightness of the results and simplicity of 
computation is present. 

E. Reduced transition matrix for infinite population 

The model presented is also adaptable to the case of infinite 
population. In fact the probability transition matrix P can be 
adapted by considering 


D. Absorbing state validity 

Since reducing the transition matrix by means of an absorb¬ 
ing state is an approximation, we now aim at verifying the 


pa = J2J2 


X t e 


-A 


t =0 6=0 


t\ 


Prtt-Pr) 


N B -b 


q(s\t + b) 


( 11 ) 

















Although in this case P has infinite size, it has been just 
shown that is possible and practical to reduce the P matrix 
by means of an absorbing state. With this expedient, also this 
case can be treaten since P is transformed to an equivalent 
matrix with finite size. 

VI. Design settings 

In this section we will resume the role of some crucial 
parameters such as M, po , p r and explain how the commu¬ 
nication state changes when changing them and what are the 
relations to be considered when designing DVB-RCS2 com¬ 
munications in RA mode. In fact, since optimization depends 
on constraints and degrees of freedom of the particular study 
case, a deep investigation of the relation among these crucial 
parameters is needed in order to understand tradeoffs that have 
to be faced. Moreover, the same discussion is introductory 
to understand how control policies can help to obtain better 
performance and a stable channel even in case of statistical 
fluctuations that could yield to instability. 


A. Degrees of freedom 


Consider Figure [10| The first thing we can notice is that 
M and po only influence the channel load line while chang¬ 
ing p r corresponds to modifying the equilibrium contour. In 
Figure |10(a)| we can see that po influences the slope of the 
channel load line that becomes steeper when increasing p 0 . 
This is intuitive to und erstand since for fixed M, if p$ < p$ 
then < G^. Figure 10(b)| shows that M shifts the channel 
load line up and down if the population size is respectively 
increasing or decreasing. From Equation [3] we can notice that 
the value M is nothing else that what is known in literature as 
q , i.e. the intersection of the line with the y-axis. Finally p r 
determines a shift upward of the channel load line for smaller 
values of p r (Figure |10(c)| ). As already shown, if the value 
p r is sufficiently small, it is possible to stabilize a channel 
initially unstable. This represents the ground base for one of 
the control policies shown in the next section. 


B. Constraints 


Concerning possible constraints , we can identify the fol¬ 
lowings 

1) stability 

2) maximize po or guarantee a value po > p™ in 

3) maximize M or guarantee a value M > M mm 

4) minimize the delay or guarantee that a certain value 
Av[D pkt ] max is not exceeded 

5) maximize the throughput or guarantee a value 

\ r^min 
t*OUT ^ D OUT 

For constraint 1) it has been previously shown that the 
stability corresponds to having a single intersection between 
the channel load line and the equilibrium curve (i.e. a single 
point of equilibrium) before the throughput peak. Constraints 
from 2) to 5) can be graphically represented on the ( GtMb ) 
plane by the following formulas: 


2 ) 


G t -N s 
M -N b 


>Po 


( 12 ) 



(a) Changing p 0 



(b) Changing M 



(c) Changing p r 


Figure 10. Graphical representation of the changes in the {Gt,Nb) plane 
when increasing or decreasing one of the parameters’ values. 


3) 


4) 


5) 


Gt-Ns 
Po 

N b 


G t -N] 

G t > G, 


+ N b > M mm 


< Av[Dpht\ 


min 

OUT 


(13) 

(14) 

(15) 


In other words, each formula simply represents an admis¬ 
sible area where the operational point has to be, as illustrated 
in Figure |TT] 
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(a) po constraint 



(b) M constraint 


the intersection of the respective admissible areas. Therefore 
we can conclude that the designed model not only allows to 
forecast stability and throughput of the communication, but 
also to easiliy represent its constraints. This enables a design 
that does not need any use of complicated formulas in order 
to understand feasible settings. 

Notice that the choice of N s has not been considered as 
a design parameter. This choice have basically two reasons: 
first of all in DVB-RCS2, the number of slots allocated for 
RA is dynamically decided by the NCC; secondly, the best 
choice for N s can always be considered between 100 and 200 
since a lower value would degrade too much the throughput 
performance while a bigger value would degrade too much the 
delay performance. 

VII. Packet degree and RA technique 

In the previous sections we have concentrated on the best 
selection of po , p r and M once the number of replicas per 
packet (packet degree) is chosen. However also the number 
of packet copies is a design parametei^] For this reason, this 
section analyzes the choice on the number of copies per packet 
to be sent. Moreover, a comparison on the advantages and 
disadvantages of CRDSA with regard to SA are outlined. 



(c) Delay constraint 



(d) Throughput constraint 


Figure 11. Graphical representation of the admissible areas for various 
constraints. 



Figure 12. Packet Loss Ratio for different packet degrees when Ns = 100 
slots 

Figure [12] shows that up to moderate channel loads, a 
greater packet degree is always convenient not only in terms 
of PLR but also in terms of throughput and packet delay. But 
considering stability to be one of our constraints, we can see 
from Figure [13] that for p r = 1 SA offers the most reliable 
performance since after the throughput peak any CRDSA 
curve rapidly degrades with N B while SA degrades in a slow¬ 
paced manner. As a result, for moderate channel load lines 
with Gt < 0.4 for N B = 0, SA is able to support up to 


5 In this paper we compare various packet degrees under the assumption 
of equal received power per replica. For discussion on the power efficiency 

Any study case will be characterized by one or more degrees when considering transponder’s average power limits refer to (T6). Moreover 
r r i i i • t r- the additional overhead due to the need of more pointers when more copies 

of freedom and one or more design constraints . In case of more are sent ^ as b een omitted, since it can be considered neglibile with regard to 

than one design constraint, the resulting constraint will be the burst size especially when few copies (no more than 4 or 5) are considered. 
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Figure 13. Equilibrium contour for different packet degrees when N$ = 100 
slots 


approximately 450 users ensuring stability, while for CRDSA 
with 2 replicas the admissible population is approximately 
50% less and for CRDSA with 3 replicas the admissible M 
is less than 150. 

However, as previously said, the operational point for 
CRDSA always occurs in a point with higher throughput, 
smaller backlogged size and smaller packet delay. Moreover 
in Section [V] it has been shown that depending on the channel 
load line the actual FET could be so big that instability could 
be accepted. For example, for M = 450 and po = 0.089 the 
computed FET is practically zero, in the sense that even though 
we know that a certain probability exists, it is so small that 
the computation resulted to be zero due to its littleness. 

Considering M = 250 and po = 0.2 SA presents a globally 
stable behaviour around its throughput peak ( Gout = 0.36) 
for Ng = 65 while CRDSA with 3 replicas presents a locally 
stable equilibrium point with throughput equal to 0.5 , TVj 
close to 0 and unstable equilibrium point for Ng « 43. 

Using formulas in Section [V] the obtained average FET 
is approximately 10 4 frames. In this case, if we are willing 
to maintain the advantages of CRDSA with 3 replicas (e.g. 
its lower PLR) some countermeasures to ensure stability 
are required. In the next Section we will show simple yet 
effective control policies able to counteract the event of drift to 
saturation. After that, the convenience of using such a dynamic 
policy in the case of finite and infinite users’ population will 
be discussed. 


VIII. Control Limit Policies 

If we are in the case in which the considered channel is 
not stable, some kind of solution to counteract the possibility 
of drift to saturation is required. Two straightforward solutions 
could be: use a smaller value for the retransmission probability 
giving then rise to a larger backlogged population for the same 
throughput value; allowing a smaller user population size M 
thus resulting in a waste of capacity. 

A third solution is the use of control limit policies to control 
unstable channels by applying the countermeasures above in a 


dynamic manner. In this section we analyze control procedures 
of the limit type ED able to ensure stability in the case of 
stationary parameters. However the same reasoning can also 
be applied in the case of dynamic ones (such as non-stationary 
M) with the necessary modifications. 

Two simple yet effective dynamic control procedures are 
considered: the input control procedure (ICP) and the re¬ 
transmission control procedure (RCP). These two control 
procedures are based upon a subclass of policies known as 
control limit policies, in which the space of the policies is 
generally composed of two actions plus a critical state that 
determines the switch between them, known as control limit. 
In this case the control limit is a critical threshold for a certain 
number of backlogged users Nb - 

A. Input Control Procedure 

This control procedure deals with new packets to transmit. 
In particular, two possible actions are possible: accept (action 
a) or deny (action d ) and the switch between them is deter¬ 
mined by the threshold Nb as previously mentioned. 

B. Retransmission Control Procedure 

As the name says, the retransmission control procedure 
deals with packets to retransmit and in particular with their 
retransmission probability. In particular two different 
retransmission probabilities p r and p c are defined 
that represent respectively the action taken in normal 
retransmission state (action r) and in critical state (action c). 
From the definition above it is straightforward that it must be 
p r > p c - The switch between these two modes is determined 
once again by the threshold Nb- 




(a) ICP (b) RCP 

Figure 14. Control limit policy examples. 


Figure |T4] graphically represents these 2 cases. As we can 
see both policies accomplish the same task of ensuring channel 
stability. However, while ICP controls the access of thinking 
users, RCP controls the access of backlogged users. 

IX. On the use of dynamic control policies for 

FINITE POPULATION 

In this section we discuss the result of the application of 
dynamic control policies for CRDSA with 2 and 3 replicas 
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compared to SA and to the case of static design (p r = p c ) 
when a finite population is assumed. Consider the case pre¬ 
viously analyzed in Section |VII| and reported in Figure [T5| 
with M = 250, po = 0.2 and p r = 1. Moreover we assume 
p c = 0.75 for CRDSA with 2 copies and p c = 0.53 for 
CRDSA with 3 copies, that are the biggest values for p c that 
ensure stability in the two respective cases. 



Figure 15. Equilibrium contour for SA and CRDSA with 2 and 3 replicas 
in normal and critical RCP state 



Figure 16. Throughput over Nb for dynamic retransmission policies 

Figures [16] and [IT] show the results in terms of throughput 
and average packet delay when varying the threshold of the 
control limit policy Nb • As we can see the results among 
the various cases are similar in terms of maximum throughput 
and minimum delay achieved, except for ICP with 3 replicas 
that shows slightly higher throughput and lower delay than the 
other cases. 

For SA instead, the resulting channel for M = 250 and 
Po =0.2 does not need any dynamic control policy since it is 
globally stable with throughput equal to 0.369 and delay equal 
to 1.76 frames. It is now clear that despite instability that needs 



Nb 


Figure 17. Packet delay over Nb for dynamic retransmission policies 


to be controlled, the time spent in critical state is so small that 
results shown for CRDSA with 2 and 3 replicas are better than 
SA. Therefore we can affirm that CRDSA always gets better 
performance than SA despite instability. 

It can also be noticed that in any case considered here, there 
is a certain range of flatness around Ng. This range changes 
depending on the policy applied, the number of replicas, the 
time spent in critical state and the performance degradation 
between normal and critical state. The width of this flat 
region is really important since inreal satellite communication 
scenarios large propagation delay has to be considered. 

A. Control policies in case of propagation delay 

Differently from the case of immediate feedback assumed so 
far throughout the paper, when the Nb threshold is crossed the 
control limit policy needs a certain propagation time from the 
NCC to the satellite terminals before it is applied. As a result 
the effective threshold is not Nb but some value Nb+8, where 
S depends on the drift of the number of backlogged users for 
the particular case and on the propagation delay. Therefore the 
wider is the flat region, the less the performance optimality of 
the controlled channel will be influenced by the propagation 
delay [17]. 

Moreover, we can see that in the case of 2 replicas a 
degradation of the throughput is also found for Nb < N%. In 
general terms, when a control limit policy is applied it does not 
make sense to choose Nb < A^, therefore this result does 
not seem to be of interest for our analysis. However, when 
Nb < Nb the control limit policy switches back to normal 
state, but as previously mentioned users do not immediately 
receive the acknowledgment for this new change due to the 
propagation delay so that the critical policy might actually 
be left for some value Nb ~ 8 resulting in the possibility of 
suboptimal performance (see ICP with 2 copies per packet). 

B. On the use of dynamic control policies 

Let us now discuss the results when instead of dynamic 
control limit policy, the channel is designed by statically 








































































decreasing the retransmission probability so that the channel is 
always stable. Differently from SA, in this case decrementing 
p r has really small influence on the overall performance 
especially when the first point of intersection between the 
equilibrium contour and the channel load line takes place 
before the throughput peak and sufficiently far from it. As 
a matter of fact with p r = 0.53 for CRDSA with 2 replicas 
and with p r = 0.75 for CRDSA with 3 replicas the following 
results are obtained: 

• CRDSA with 2 replicas: Gout = 0.49; Av[D pkt \ = 1.01 

• CRDSA with 3 replicas: Gout = 0.5; Av[D p kt] = 1 

In conclusion we can say that with stationary values of M 
and po , in most of the cases the use of dynamic policies 
for CRDSA does not appear to be a convenient solution 
compared to a static design if a choice can be made. This 
can be affirmed also in light of the more realistic case of non- 
immediate feedback in which a certain performance degrada¬ 
tion is present. The only case in which dynamic policies could 
be really necessary with stationary M and po is when M is 
so big that the required p r rendering the channel stable really 
degrades the performance thus requiring the use of dynamic 
policies. Nevertheless the analysis above is useful for future 
work on the case of non-stationary po and M that does not 
allow to a-priori set the retransmission probability in order 
to be sure that the channel is always stable. Moreover the 
analysis presented for the case of finite M is introductory to 
understand the application of dynamic policies to the case of 
infinite population. 

X. On the use of dynamic control policies for 

INFINITE POPULATION 

Finally we discuss the use of dynamic control policies in 
order to provide stability in the case of infinite population. 
As a matter of fact this case is closer to the real DVB-over- 
satellite application scenario in which a big number of users 
is expected to generate new traffic independently of whether 
previous transmissions are still pending or not, so that actual 
channel traffic can be modeled as a Poisson process. 

In the case of infinite population, ICP is the only policy that 
can ensure stability since decreasing p r could not be sufficient 
to recover stability. In fact, RCP might be used as a mean to 
recover the channel communication from instability by setting 
p c so that the operational point falls again in the region of 
stability. However such a procedure does not always ensure 
stability from a theoretical point of view since due to the 
satellite propagation delay, the switch from and to critical 
state need some time to propagate meanwhile the number 
of backlogged users will continue to rapidly and indefinitely 
increase so that at the end the chosen value for p c could not 
be anymore sufficient. 

Figure |~j~8j and |T9| show throughput results depending on the 
number of new packet arrivals A at each frame. As we can 
see the bigger the A value, the higher the throughput but also 
the narrower the flat region. The reason for such a result is 
similar to what has been explained in the previous sections for 
the finite population case. In particular we can see from Table 
II that the closer the operational point is to the throughput 
peak, the smaller is the value of Ng. 



Figure 18. Throughput over Nb for ICP in case of infinite population 



Figure 19. Packet delay over Nb for ICP in case of infinite population 



Figure 20. Percentage of time spent in critical state over Nb for ICP 
case of infinite population 
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K = n 

Ng = 18 

NU = 25 

NU = 32 


Table II 

Points of unstable equilibrium for CRDSA with 3 copies 
DEPENDING ON A 


Figure [19] shows that obtained results concerning delay 
reflect throughput results since the wider is the flat region 
for the throughput, the wider is the flat region with delay 
equal to 1. Finally Figure [20] shows that for bigger values 
of A users spend some time in critical state also in the flat 
throughput region giving place on one hand to the possibility 
of suboptimality due to the propagation delay and on the other 
hand increasing the overall PLR by reason of rejected packets. 

XI. Conclusions 

In this paper a model for design of a congestion-free 
system using Contention Resolution Diversity Slotted Aloha as 
Random Access mode in DVB-RCS2 has been outlined. This 
model allows to study and predict stability and to calculate 
the expected throughput as well as the expected delay at the 
channel operating point. Moreover, this kind of representation 
allows a simple yet effective view of how the communication 
behaves when changing one or more of the key parameters 
thus allowing a design that can easily consider the tradeoff 
among multiple constraints. 

A simplified framework for the calculation of the First Exit 
Time has also been introduced and its importance for the 
evaluation of unstable channels has been highlighted. All the 
presented models and tools can be extended to consider a 
different frame size, power unbalance, the introduction of First 
Error Correction or any other modification. The central part 
of the paper discusses design settings and the convenience of 
using CRDSA over SA. 

In the last part, dynamic retransmission policies of the 
control limit type have been introduced and their effectiveness 
and convenience have been discussed against static design. 
In particular, concerning the case of finite user population, 
found results have shown that differently from SA in which 
dynamic retransmission policies are generally preferred, in 
the case of CRDSA a static design that ensures stability (i.e. 
by decreasing the retransmission probability for backlogged 
packets) is generally more convenient. Finally, dynamic re¬ 
transmission policies have been applied to the case of infinite 
user population. Found results have shown that control limit 
policies are able to stabilize the channel also in the case of 
infinite user population while ensuring a performance close 
to the best achievable even in presence of large propagation 
delay as it is the case in a geostationary satellite system used 
for DVB-RCS communications. 

References 

[1] J. Maisonneuve, M. Deschane, J. Heiles, Wei Li, Hong Liu, R. Sharpe, 
and Yiyan Wu, ”An Overview of IPTV Standards Development,” IEEE 
Trans. Broadcast., vol. 55, no. 2, pp.315-328, Jun. 2009. 


[2] Digital Video Broadcasting (DVB): Second Generation DVB Interactive 
Satellite System (DVB-RCS2); Part 2: Lower Layers for Satellite standard, 
ETSI standard EN 301 545-2 VI.1.1 , Jan. 2012. 

[3] L.G. Roberts, ’’ALOHA packet systems with and without slots and 
capture”, ARPANET System Note 8 (NIC11290), Jun. 1972. 

[4] N. Abramson, ’’The throughput of packet broadcasting channels”, IEEE 
Trans. Commun., vol.25, pp. 117-128, Jan. 1977. 

[5] E. Casini, R. De Gaudenzi, and O. del Rio Herrero, ’’Contention Resolu¬ 
tion Diversity Slotted ALOHA (CRDSA): An Enhanced Random Access 
Scheme for Satellite Access Packet Networks”, IEEE Trans. Wireless 
Commun., vol.6, no.4, pp. 1408-1419, Apr. 2007. 

[6] G.L. Choudhury and S.S. Rappaport, ’’Diversity ALOHA - A random 
access scheme for satellite communications”, textitlEEE Trans. Commun., 
vol.31, pp.450-457, Mar. 1983. 

[7] M.J. Hagh and M.R. Soleymani, ’’Application of Raptor Coding With 
Power Adaptation to DVB Multiple Access Channels,” IEEE Trans. 
Broadcast., vol.58, no.3, pp.379-389, Sep. 2012. 

[8] M.J. Hagh and M.R. Soleymani, ’’Constellation Rotation for DVB Multi¬ 
ple Access Channels With Raptor Coding,” IEEE Trans. Commun., vol.59, 
no.2, pp.290-297, Jun. 2013. 

[9] A. Meloni and M. Murroni, ’’CRDSA, CRDSA++ and IRSA: Stability and 
performance evaluation”, in Proc. 6th ASMS and 12th SPSC Workshop, 
Baiona, Spain, Sep. 2012, pp. 220-225. 

[10] L. Kleinrock and S. Lam, ’’Packet Switching in a Multiaccess Broadcast 
Channel: Performance Evaluation”, IEEE Trans. Commun., vol.23, no.4, 
pp.410-423, Apr. 1975. 

[11] S. Lam and L. Kleinrock, ’’Packet Switching in a Multiaccess Broad¬ 
cast Channel: Dynamic Control Procedures,” Communications, IEEE 
Trans.Comm. , vol.23, no.9, pp.891,904, Sep 1975 

[12] G. Liva, ’’Graph-Based Analysis and Optimization of Contention Reso¬ 
lution Diversity Slotted ALOHA”, IEEE Trans. Commun., vol.59, no.2, 
pp.477-487, Feb. 2011. 

[13] A. Ijaz, A.B. Awoseyila, and B.G. Evans, ”Signal-to-Noise Ratio Es¬ 
timation Algorithm for Advanced DVB-RCS Systems”, IEEE Trans. 
Broadcast, vol.58, no.4, pp.603-608, Dec. 2012. 

[14] Dae-Ki Hong and Sung-jin Kang, ’’Joint frequency offset and carrier 
phase estimation for the return channel for digital video broadcasting,” 
IEEE Trans. Broadcast , vol.51, no.4, pp.543-550, Dec. 2005. 

[15] L. Kleinrock, ’’Queueing Systems: Volume I Theory; Queueing Sys¬ 
tems: Volume II Computer Applications”, Wiley Interscience, New York, 
1975-1976. 

[16] A. Meloni and M. Murroni, ’’Average Power Limitations in Sliding 
Window Contention Resolution Diversity Slotted Aloha”, in Proc. 9th 
International Wireless Communications and Mobile Computing Confer¬ 
ence, Cagliari, Italy, Jul. 2013, pp. 7-12. 

[17] A. Meloni and M. Murroni, ’’Random Access Congestion Control 
in DVB-RCS2 Interactive Satellite Terminals”, in Proc. IEEE Interna¬ 
tional Symposium on Broadband Multimedia Systems and Broadcasting, 
Uxbridge, West London, England, Jun. 2013, pp. 1-6. 


Alessio Meloni received his M.Sc. degree in elec¬ 
tronic engineering from the University of Cagliari 
(Italy) in December 2010, after spending one year 
at the Czech Technical University (Prague) during 
which he developed his thesis. Since March 2011, 
he has been PhD student at the Multimedia and 
Communications Lab of the Department of Electri¬ 
cal and Electronic Engineering (DIEE) - University 
of Cagliari. From April 2011 to October 2011 he has 
been guest PhD student at the German Aerospace 
Center, Oberpfaffenhofen, Germany. His research in¬ 
terests include wireless communications, advanced random access techniques 
and satellite links. He is a member of IEEE, IEE ComSoc, CNIT and GOLD 
member of BTS. 







14 


Maurizio Murroni (M’Ol) graduated (M.Sc.) with 
honors (Summa cum Laude) in Electronic Engi¬ 
neering in 1998 at the University of Cagliari and 
in the same year he received an award for his 
thesis from Telecom Italia, inc. He has been an 
Erasmus visiting student at CVSSP Group (Prof. 
Maria Petrou), School of Electronic Engineering, In¬ 
formation Technology and Mathematics, University 
of Surrey, Guildford, U.K. in 1998 and a visiting 
Ph.D. student at the Image Processing Group (Prof. 
Yao Wang), Polytechnic University, Brooklyn, NY, 
USA, in 2000. In 2001 he received his PhD degree in Electronic Engineer¬ 
ing and Computers, from the University of Cagliari. In 2002 he became 
Assistant Professor of Communication at the Department of Electrical and 
Electronic Engineering (DIEE) of the University of Cagliari. Since the 1998, 
he contributes to the research and teaching activities of the Multimedia 
Communication Lab (MCLab) at DIEE. In 2006 he has been visiting professor 
at the Dept, of Electronics and Computers at the Transilvania University of 
Brasov in Romania and in 2011 at the Dept. Electronics and Telecommu¬ 
nications, Bilbao Faculty of Engineering, University of the Basque Country 
(UPV/EHU) in Spain. Since October 2010 he is coordinator of the research 
unit of the Italian University Consortium for Telecommunications (CNIT) 
at the University of Cagliari. His current research focuses on Cognitive 
Radio system, signal processing for radio communications, multimedia data 
transmission and processing. Dr. Murroni is a member of IEEE, IEEE Com 
Soc, IEEE BTS, IEEE DySPAN-SC and 1900.6 WG. 




